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Abstract. We present a linear analysis of the vertical structure and growth of 
the magnetorotational instability in weakly ionised, stratified accretion discs. The 
method includes the effects of the magnetic coupling, the conductivity regime of 
the fluid and the strength of the magnetic field, which is initially vertical. The 
conductivity is treated as a tensor and assumed constant with height. 

The Hall effect causes the perturbations to grow faster and act over a much 
more extended section of the disc when the magnetic coupling is low. As a result, 
significant accretion can occur closer to the midplane, despite the weak magnetic 
coupling, because of the high column density of the fluid. This is an interesting 
alternative to the commonly held view that accretion is relevant mainly in the surface 
regions of discs, which have a better coupling, but a much lower fluid density. 

Keywords: accretion discs, instabilities, magnetohydrodynamics 



1. Introduction 

The magnetorotational instability (MRI) (Balbus k Hawley 1991(BH91 
hereafter), Hawley & Balbus 1991) transports angular momentum radi- 
ally outwards in accretion discs through the distortion of the magnetic 
field lines that connect fluid elements. In protostellar discs, low conduc- 
tivity is important, specially in the inner regions (Gammie 1996; Wardle 
1997). As a result, low k modes are relevant and vertical stratification 
is a key factor of the analysis. However, most models of the MRI in 
these environments have adopted either the ambipolar diffusion or resis- 
tive approximations and have not simultaneously treated stratification 
and Hall conductivity. We present here a linear analysis of the MRI, 
including the Hall effect, in a stratified disc. 



2. Formulation 



We write the equations of non-ideal MHD about a local keplerian frame 
corotating with the disc at the angular frequency 17. The velocity field 
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is expressed as a departure from exact keplerian motion, v = V — mk , 
where V is the velocity in the standard laboratory frame (r, (p, z) and 
vk is the keplerian velocity at the radius r. We assume the fluid is 
weakly ionised, meaning that the abundances of charged species are 
low enough that their inertia and thermal pressure can be neglected, 
together with the effects of ionisation and recombination processes in 
the neutral gas. Under these conditions, separate equations of motion 
for the charged species are not required and the fluid equations are: 

^ + V-(pv)=0, (1) 



V c J X !B 

— + (v • V)v - 2nv^r + ^nvrfp --^r+ ^Vp + V$ = , (2) 

ot r p cp 

_ = Vx(vxB)-cVxE'-|OB^^. (3) 

In the equation of motion (2), $ is the gravitational potential for a 
non self-gravitating disc and ^'x/r is the centripetal term generated 
by keplerian motion. Coriolis terms 2r2i;^f and ^Ouri^ are associated 
with the use of a local keplerian frame and Cg is the isothermal sound 
speed. In the induction equation (3), E' is the electric field in the frame 
comoving with the neutrals and |f2Br^ accounts for the generation of 
toroidal by differential rotation. Additionally, the magnetic field must 
satisfy the constraint V • B = and the current density must satisfy 
Ampere and Ohm's laws. Following Wardle & Ng (1999), Wardle 1999 
(W99) and references therein, the current density is expressed as, 

J = (7 • E' = (7||Ej| +criBxE'j_ -|-(72E'_i_ , (4) 

where Ej| and E'j^ are the components of E' parallel and perpendic- 
ular to B. Note that the conductivity is a tensor with components (T||, 
parallel to the magnetic field, ui, the Hall conductivity, and C72, the 
Pedersen conductivity. 

Our model includes vertical stratification, but it is local in the radial 
direction. The vertical density distribution in hydrostatic equilibrium 
is p/po = exp(— 2:^/2iJ^) where po is the midplane gas density and H = 
Cg/^ is the scaleheight of the disc. Equations (1) to (4) were linearised 
about an initial steady state where J = v = E' = and B = Bz. Note 
that as E' vanishes in the initial state, the changes in the conductivity 
due to the MRI are not relevant in this linear formulation and only 
the unperturbed values of a are required. Taking perturbations of the 
form q = qo + 5ci{z)e^'^^ and assuming k = kz we obtained a system 
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of ordinary differential equations (ODE) in (5E (the perturbations of 
the electric field in the laboratory frame), 5B, and the growth rate 
V = iu/Q. Three parameters control the evolution of the fluid. (1) 
va/cs, the ratio of the Alfven speed and the isothermal sound speed 
of the gas at the midplane, which is a measure of the strength of the 
magnetic field. (2) Xo, the ratio of the frequency above which flux- 
freezing breaks down and the dynamical frequency of the disc. When 
Xo < 1 the magnetic field is poorly coupled to the disc. (3) ai/a2, 
the ratio of the conductivity terms perpendicular to the magnetic field, 
which characterises the conductivity regime of the fluid. 

We present here results associated with the ambipolar diffusion limit 
((Ji = 0), both Hall limits ((72 = 0, aiB^ > and aiB^ < 0) and the 
full conductivity regime (o"i = cj2)- Ambipolar diffusion is dominant 
at relatively low densities, when the magnetic field is frozen into the 
ionised component of the fluid. Conversely, Hall limits are predominant 
at intermediate densities, and are characterised by a varying degree of 
magnetic coupling amongst charged species. Typically electrons are well 
coupled, but ions and especially grains tend to be tied to the neutrals 
via collisions. We note that the Ohmic regime, which dominates at low 
densities, when the charged species are attached to the neutrals via 
collisions, is identical to the ambipolar diffusion regime for the modes 
of interest here. We integrated this system of equations vertically as 
a two-point boundary value problem for coupled ODE with boundary 
conditions SBr = SB^ = and SE'^. = 1 at z = and 5Br = SB^ = 
at z/H = 5. We refer the reader to Salmeron & Wardle (in press), for 
full details of the method and results. 



3. RESULTS 

We compare the structure and growth rate of MRI perturbations as a 
function of the magnetic coupling Xo for different conductivity regimes 
and va/cs = 0.1 (Fig. 1). At very high magnetic coupling {xo ~ 100), 
ideal MHD holds and all conductivity regimes are alike (see leftmost 
column of the figure). For a weaker coupling, such that Xo > va/cs, the 
structure of the ambipolar diffusion limit and full conductivity regimes 
is similar, signalling that ambipolar diffusion is dominant in this region 
of parameter space. Note that both perturbations peak at the node 
closest to the surface. This occurs because the maximum growth rate 
of ambipolar diffusion perturbations increases with the local x (W99), 
which in turn is a function of height. As a result, at higher z the local 
growth of the instability increases, driving the amplitude of these global 
perturbations to increase. Hall perturbations, on the other hand, peak 
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at the node closest to the midplane, as i^max is the same for all x 
(W99) in this case and the instability is not driven from any particular 
vertical location. Finally, for Xo < va/cs (rightmost column), ambipolar 
diffusion is damped (W99) and the MRI is driven by the Hall effect. In 
this case ambipolar diffusion and full conductivity perturbations differ, 
the structure of the latter ones resembling the Hall limit, as expected. 

We also explored the growth rate of the perturbations in parameter 
space. Fig. 2 (left panel) shows the growth rate of the most unstable 
perturbations {vmax) as a function of Xo for ambipolar diffusion and 
both Hall limits with vj^jcs = 0.1. At good coupling the instability 
grows at ~ 0.75O in all cases. As Xo diminishes, I'max is reduced at a 
rate that depends on the conductivity regime. The maximum growth 
rate departs significantly from the ideal value for Xo ~ 0.1 (ambipolar 
diffusion regime) and for Xo ~ 0.01 (Hall (JiBz > limit). The Hall 
aiBz < case could not be analysed for Xo < 2 because in this region 
of parameter space all wavenumbers grow (W99) and our code fails 
to converge. These results arc consistent with the finding that am- 
bipolar diffusion modes grow when x^va/cs while Hall cases require 
X > va'^/cs'^ (W99). As a result, for Xo < va/cs ambipolar diffusion 
perturbations have negligible growth, but when Hall condTictivity is 
present the instability still grows at = 0.2 — 0.3 (see Fig. 1). In fact, 
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Figure 2. Mciximum growth of the MRI. Left: as a function of Xo with va/cs = 0.1. 
Right: as a function of va/Cs for Xo = 10 and Xo = 2 (Hall UiS^ < regime only) 



we found that when Hall terms dominate, u ^ 0.75 for Xo ~ 10~^ even 
at va/cs ~ 0.01, as expected. 

Fig. 2 (right panel) shows the maximum growth rate as a function 
of va/cs for Xo = 10 and 2 (the latter one for Hall aiBz < regime 
only). For very weak va/cs, the MRI grows at close to the ideal rate 
in all regimes. Generally, increasing the strength of the magnetic field 
has little effect on Vmax until va/cs ~ 1, when all perturbations are 
damped (e.g. see Fig. 2, Xo = 10). At this va/cs, the wavelength of the 
most unstable modes becomes ~ H, the scalcheight of the disc (BH91). 
An exception to this occurs in the Hall aiBz < regime for Xo = 2. In 
this case we found unstable modes for va/cs up to 2.9. We know from 
the local analysis (W99) that once the local x < 2, unstable modes 
exist for every Uva/^ in this regime. As a result, even for suprathcrmal 
fields {va/cs > 1), there are still unstable modes with kH < 1 growing 
within the disc. 



4. Discussion 



In a real disc, different conductivity regimes are expected to dominate 
at different heights. Despite this, it is common to neglect the contribu- 
tion of Hall terms in studies of non-ideal accretion discs, which generally 
adopt either the ambipolar diffusion or resistive approximations. To 
illustrate the way Hall conductivity can affect the dynamics and evolu- 
tion of the fluid we modelled the MRI with a set of parameters such that 
the Hall regime dominates closer to the midplane while ambipolar diffu- 
sion is predominant near the surface, as expected to occur in a real disc. 
We compared these results with the ambipolar diffusion approximation 
(see Fig. 3). Note that in the full conductivity case the extent of the 
dead zone is reduced and the growth rate is increased in relation to the 
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Figure 3. Structure and growth rate of the MRI when the Hall regime is dominant 
close to the midplanc and ambipolar diffusion dominates near the surface (left panel) 
and under the ambipolar diffusion approximation (right panel), va/cs = 0.01. 

ambipolar diffusion limit. Tfiis fias implications for the global evolution 
of the disc, specially considering that even though the MRI is damped 
in the dead zone, velocity fluctuations persist, driven by turbulence 
in the active zones, and transport angular momentum through non- 
axisymmctric density waves (Stone &: Fleming 2003). On the other 
hand, a quiescent zone would allow dust grains to settle towards the 
midplane and begin to assemble into planetesimals (e.g Weidenschilling 
& Cuzzi 1993), so the process of planet formation could also be affected 
by Hall conductivity. 

Clearly, more detailed modelling is required. The assumption of con- 
stant conductivity is unrealistic and will affect the nature of the global 
modes. An analysis including a z-depcndcnt conductivity is underway 
to examine more fully the MRI in low conductivity discs and to quantify 
the importance of Hall terms. 
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